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Abstract: The Jaynes-Cummings-Gaudin model describes a collection of n spins coupled to an 
harmonic oscillator. It is known to be integrable, so one can define a moment map which associates 
to each point in phase-space the list of values of the n + 1 conserved Hamiltonians. We identify all 
the critical points of this map and we compute the corresponding quadratic normal forms, using the 
Lax matrix representation of the model. The normal coordinates are constructed by a procedure 
which appears as a classical version of the Bethe Ansatz used to solve the quantum model. We 
show that only elliptic or focus-focus singularities are present in this model, which provides an 
interesting example of a symplectic toric action with singularities. To explore these, we study in 
detail the degeneracies of the spectral curves for the n = 1 and n = 2 cases. This gives a complete 
picture for the image of the momentum map (IMM) and the associated bifurcation diagram. For 
n = 2 we found in particular some lines of rank 1 which lie, for one part, on the boundary of the 
IMM, where they behave like an edge separating two faces, and which go, for another part, inside 
the IMM. 

1 Introduction. 

Suppose we have a Hamiltonian system with Hamiltonian H on a phase space of dimension 2n with 
coordinates x\, ■ ■ ■ ,X2 n - A critical point x^ is an equilibrium point 



dH 

dxj 



= 0, j = l,---,2n 

_ (o) 

Xj —Xj 



We can expand the Hamiltonian around such a point, Xj = Xj + 5xj, getting a quadratic form 

H = H<®+Q + ---, Q = ^2Q ij 5x i 5x j 

The question then arises to simplify or "diagonalize" the quadratic form Q. This problem is 
simple if we allow Sl2 n (R) transformations of the 5xj, the matrix Qij can then be diagonalized in 
the orthogonal group. However in the present case the natural group is the group of symplectic 
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transformations and the diagonalization problem is not straightforward. Its solution has been 
given by Williamson [H [2] in terms of six families of elementary quadratic forms on which we can 
decompose H. Each family corresponds to a certain pattern of eigenvalues of the associated linear 
Hamiltonian flow. There are two families for the eigenvalue 0, a third one for a pair of opposite real 
eigenvalues (hyperbolic flow), two more for a pair of purely imaginary eigenvalues (elliptic flow), 
and the last one for quartets of complex eigenvalues of the form ±a =t i/3. Besides diagonal blocks, 
each family admits non-diagonal Jordan blocks of unlimited size. 

The above problem of reducing H to its normal form is not completely trivial. Even less trivial 
is the case of an integrable system, in which case we have n Poisson commuting Hamiltonians Hi. 
These functions define the moment map, from R 2n to R n , which sends x to (H\(x), H n (x)). A 
critical point is a simultaneous critical point for all the Hi, in other words, it is a point at which 
the differential of the moment map vanishes. Expanding the Hi around such a point, we get n 
Poisson commuting quadratic forms Qi. For each of them, Williamson theorem applies, but since 
they commute, we can reduce them to normal forms simultaneously. The simplest case is of a 
nondegenerate singularity, for which the Hamiltonian flows associated to the n quadratic forms Qi 
realize a Cartan subalgebra of the Lie algebra sp n (R). The fact that each Qi Poisson commutes 
with the others yields strong constraints on the form of the elementary blocks which can appear in 
the Williamson reduction of Qi. It turns out that integrability excludes non-diagonal Jordan blocks, 
therefore there exist canonical coordinates q\, ■ ■ ■ q n , pi, ■ ■ ■ ,p n such that the above quadratic forms 
can be reduced into the following quadratic polynomials: [3j 

^elliptic = p 2 + g 2 ; j = 1)2 ,..., mi 
^hyperbolic • . 

P- JP = pjqj, j = mi + 1, • • • mi + m 2 

pfocus-focus = p . q . +p . +iq . +li j=mi+m2 + 1 ... 
pfocus-focus = p . q . +1 _ p . +iq . } ... mi+m2 + 2m3 

where mi + m<i + 2m-3 = n. The triple (mi, m-i^m^) is called the type of the critical point. It 
encodes all the qualitative information on the system, in particular on the dimensions of fibers of 
the moment map in the vicinity of the corresponding critical value. The above set of quadratic 
forms is refered to as the normal form for a non-degenerate critical point of an integrable system. 
In the focus-focus case, pj' ocus - focus anc i pfocus-focus can bg com bi nec i i n t Q a single complex 

quantity. Setting 

B i = Pj + Wj+ii B j+i = Qj + 



we have 



r> d rjfocus— focus • r>focus— focus 



In this paper we observe that the solution to this problem of the simultaneous reduction of the 
conserved Hamiltonians in the vicinity of critical points is the classical analog of Bethe Ansatz which 
is used to diagonalize simultaneously the quantum commuting Hamiltonians Hi. We demonstrate it 
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on the particular example of the Jaynes-Cummings-Gaudin model [U [6j [7] , but the construction 
is clearly a general one. As a result we will show that, in this model, singularities are only of two 
types: elliptic and focus-focus. 

The above critical points are the points where the moment map 

R 2n^ R n. x ->(Hi(x),...,H n (x)) 

is of rank zero. Other manifolds where the rank of the moment map is not maximal are also of great 
interest. Their images in R n by the moment map constitute the bifurcation diagram. As explained 
in particular by Michele Audin [Hj, a powerful tool to construct this diagram for an integrable model 
is the correspondence between singularities of the moment map and degeneracies of the associated 
spectral curve. We will show explicitely the results of this method for the Jaynes-Cummings-Gaudin 
model with two and three degrees of freedom. 

Another advantage of this model is that it admits one free parameter for each degree of freedom. 
The types of the critical points therefore depend on these parameters. In particular, it is possible 
to choose them in such a way that all the critical points are of the elliptic type. In this case, 
the Hamiltonian flows of the conserved Hi can be used to define a torus action on phase space. 
A famous theorem states that, if we replace the Hi's by the corresponding action variables, the 
image of the moment map is a convex polytope [H[T0]. However, for most values of the parameters, 
focus-focus singularities appear which provides examples of the more general "almost toric" action 
studied in and the important phenomenon of "monodromy" discovered in |12| . preventing the 
existence of global action- angle variables, comes into play. 

While the one-spin Jaynes-Cummings model was already re-discovered by mathematicians |13j . 
it seems that the results concerning the bifurcation diagrams of the two-spin Jaynes-Cummings 
model are new. The geometric richness of these diagrams undoubtedly calls for further studies. 

The paper is organized as follows. In section 2 we present the Jaynes-Cummings-Gaudin model 
and recall some basic facts about its Lax representation. In section 3 we study the critical points of 
the model and their normal forms are computed in section 4 using a classical analog of Bethe Ansatz. 
In section 5 we discuss the monodromy. Section 6 is a short introduction to separated variables and 
explains that the different strata of the bifurcation diagram correspond to the degeneracies of the 
spectral curve. In sections 7 and 8 we study the bifurcation diagram for the one-spin and two-spin 
models respectively. 

2 The classical Jaynes-Cummings-Gaudin model. 

This model, where a collection of n spins is coupled to a single harmonic oscillator, has been used 
for more than fifty years in atomic physics to describe the interaction of an ensemble of atoms 
with a mode of the quantized electromagnetic field [H [61 [2] . It derives from the following 
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Hamiltonian: 

n n 

H = ^(2 Cj + u)sl +oobb + ^2 (bsj + 6s+) (1) 

3=1 3=1 

The Sj are spin variables, and b, b is a harmonic oscillator. The Poisson brackets read 

{s^,s b j } = -e abc s c j , {b,b} = i (2) 
The Sj brackets are degenerate. We fix the value of the Casimir functions 
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Phase space has dimension 2(n + 1). In the Hamiltonian we have used 



which have Poisson brackets 



sf = s) ± is] 



{s*,sf} = ±isf, {s+s-} = 2is) 



J ' 3 J 3 ' 1 J ' J J J 

The equations of motion read 

dH TO ^ 

b = -i^ = -iub-iY,sJ (3) 

3=1 

.dH , .dH _ .,- _ , , , 

sj = i^s+-2i^s* = i{2e j + u)s+-2ibs* (5) 
ij = -i^j + JM^J = -<(&,- + W ) a j + 2i^ (6) 
This is an integrable system. To see it we introduce the Lax matrices 

n _ _ 

L(A) = 2A<r z + 2(kr+ + + X] 7^" (7) 

3=1 ~ ^ 

M(A) = -i\a z -i^a z -i(ba + + ba-) (8) 
where cr a are the Pauli matrices. 

i\ v fo -A 2 /i o 

1 o ' a = U ' * = -1 
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and we have defined 



a ± = ^(a x ±ia y ), [a z , a±] = ±2a ± , [a + ,a-]=a z 



It is not difficult to check that the equations of motion are equivalent to the Lax equation 

L(X) = [M(X),L(X)} 

Let 

L(X) = 



A(X) B(X) \ 
C(X) -A(X)J 



we have 

A(X) 
B(X) 
C(X) 

One has the simple Poisson brackets 

{A(X),A^)} 
{B(X),B^)} 
{C(A),C( M )} 
{A(X),B(jj,)} 

{A{\),CM} 

{B(X), C( M )} 



2A + E 



j'=i 

n 



X — 6j 



3 



= 26 + E 



j'=i 



A — e,- 









X — fi 



(B(\)-B(jjl)) 
-(C(A)-C(m)) 



A — /U 

(A(A)-A(m)) 



X — fi 



One can rewrite these equations in the usual classical r-matrix form 

fl2 



{Li(A),L 2 (/x)} 



X — fi 



,L 1 (X) + L 2 ( Ji ) 



where 



Pi 



12 



/l 0\ 

10 

10 

\0 1/ 



(9) 



(10) 

(11) 

(12) 

(13) 
(14) 
(15) 

(16) 
(17) 
(18) 
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It follows immediately that TV (L 2 (A)) = 2A 2 (A) + 2B(A)C(A) Poisson commute for different values 
of the spectral parameter: 

{Tr(L 2 (A 1 )),Tr(L 2 (A 2 ))} = 
Hence A(A) = ^Tr(L 2 (A)) generates Poisson commuting quantities. One has 



where the (n + 1) Hamiltonians i/j, j = 1, • • • , n + 1 read 

tf n+1 = bb+Yl s j ( 2 °) 
3 

and 

H j = 2e j s3 + {bs+ + bsJ) + Y,^ lL , 3 = !,■■•," (21) 

One can easily verify that, indeed, {Hi,Hj} = for i, j = 1, • • • , n+1, hence the system is integrable. 
The Hamiltonian eq.JT]) is 

n 
3=1 

3 Critical points 

The critical points are equilibrium points for all the Hamiltonians Hj, j = 1, • • • , n + 1. At such 
points the derivatives with respect of all coordinates on phase space vanish. In particular, since 

dHn+i _ , <9#j _ _ 
db ~ ' 56 ~ Sj ' 

we see that the critical points must be located at 

b = b = 0, sf = 0, s| = ejs, ej = ±1 (22) 

When we expand around a configuration eq. ([22j) . all the quantities (b, b, s~j~, sj) are first order, 
but Sj is second order because 



-2 



s 3 = e i V s '" s i s i = se i ~ t a 8 i a i + " ' ' e i = ±X 

It is then simple to see that all first order terms in the expansions of the Hamiltonians Hj vanish. 
Hence we have found 2 n critical points. 
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4 Normal Forms. 

We want to expand the Hamiltonians Hj around the equilibrium points eq. (|22p and write them in 
normal form. Symbolically : 

Hj = ^2 Ej, a a a a a (23) 

a 

where CLq> f Q>ct 9>rG independent harmonic oscillators (in the elliptic case). Note that if we quantize 
the above Hamiltonians in this approximation, their diagonalization is immediate: the normal 
coordinates a a become spectrum generating operators. Their construction must therefore be very 
much related to the simultaneous diagonalization of the Hj . But the tool to solve this problem is 
well known: Bethe Ansatz. 



Inspired by this remark, we return to eqs. (fTU1 -[T2" j) and eqs. p^1 - [TS|) . Now, we have 

j^Tr L 2 (A),CV)| = J±-(a{\)C{») - A^)C(X)) 

When we expand around a critical point, the Hamiltonians are quadratic. Remark that C(/i) is 
first order and therefore the Poisson bracket in left hand side is linear. Now A(X) is constant plus 
second order, so that in the right-hand side we can replace A(X) and A(/j,) by their zeroth order 
expression : 

n 

A(X) ~ a(X) = 2A + £ 

j = l J 

and we arrive at 

j^Tr L 2 (A),C(/,)| = ^-(a(A)C(M)-a(M)C7(A)) (24) 

Going back to eq. (|23]) . we see that the a a are such that 

{Hj,a a } = E jyCt a a (25) 

Eq. (j24p will be precisely of the form of eq. ([25p if we can kill the unwanted term C(A). This is 
achieved by imposing the condition 

a(n) = 0, "Classical Bethe Equation" (26) 

This is an equation of degree n + 1 for fi. Calling fj,{ its solutions, we construct in this way n + 1 
variables C(fii). Remark that by eq. (|15p . they all commute 

{C(^),C( H )} = (27) 



7 



Since phase space has dimension 2(n + 1) this is half what we need. To construct the conjugate 
variables, we consider eq. (|18p . In our linear approximation it reads 

2i 

{B(pi),C(jij)} = — - — (o(/ii) -a(fij)) 
Hi 

If ^ and fj,j are different solutions of eq. (f26|) , then obviously 

{B( lM ),C( H )} = 0, fa^nj (28) 

If however fj,j = [i% then 

{B{^),C{Hi)} = 2ia'{^) (29) 

Finally, by eq. (|14p we have 



{B(/ii),B(Mi)} = (30) 
Up to normalisation, we have indeed constructed canonical coordinates ! 

It is simple to express the quadratic Hamiltonians in theses coordinates: 

\tt l' W = AX) + E al( ^_ H) B( m )C M (31) 

This has the correct analytical properties in A and together with the Poisson brackets eqs. (|27!28l 
I2UI l3"Uj) we reproduce eq. (f!Mj) . Note that there is no pole at A = fj,j because a(fij) = 0. Expanding 
around A = oo we get 

H n+1 = s e k + E 2a'(fi-) B ^ C ^ 
and computing the residue at A = e,-, we find 



Hj = sej 



^ + E 



se k 



^-T 1 2a' (m) ej - Hi 



We can invert these formulae: devide eq. ([3Tj) by A — fij an d take the residue at A = fij- Since 
a(fij) = we get 

l -TrL\H ] ) = B(H 3 )C{H j ) 

or explicitly 

n 2H k " s 2 



We can now make contact with the Williamson classification theorem. 
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If [ij is real we have B(nj) = C(fj,j) and we set 

C(Vj) = \J\a'(fJ,j)\(pj + ie, qj), B(pj) = yV(/^)IOj - Qj ) 
where ej = — sign(a'(^)) and pj,qj are canonical coordinates. Then eq. (|29p is satisfied. Moreover 

B{ H )C{^) = \a' {^)\{p) + qj) 
i.e. we have an elliptic singularity. 



If jUj is complex, there is another root fij+i = /Uj which is its complex conjugate. Then B(/j,j) = 
C(//j+i). We introduce canonical coordinates pj,qj,pj + \,qj + \ and set 

-^i) = -ia'i^iQj + %+l), = Pi + ^i+i 

C^i) = Pi - C^i+i) = ia '(H+i)(Qj ~ kj+i) 

then eas. (l27l28l29l3Up are satisfied and 

B(^)C(Mi) = -ia\^)(Pf ocus ' fo ^ s + iPf^ us ~ focus ) 
B(p j+ i)CQM j+1 ) = m'( Mj+1 )(P/ ocus - /ocus -iP/;r s - /ocus ) 
i.e. we have a focus-focus singularity. 

It remains to see when the classical Bethe roots are real and when they are complex. Let us 
assume all the spins are down: tj = —1, j = 0, ■ • ■ , n — 1. The condition a(/u) = reads 

"-5?i^ (32) 

The graph of the curves y = ji and y = I are presented in Fig. [I]. In that case we have 

J r*' J 

n + l real roots. The singularity is elliptic meaning that this critical point is locally stable. This is 
the situation that prevails when ^ a < 0. 

Suppose now that e,- = 1, j = 0, ■ ■ • , n— 1. The graph of the curves y = u and y = — § — — 
are presented in Fig [2]. The situation is more complex, we can have n—X real roots and a pair of 
complex conjugated roots, or n + 1 real roots, depending on the values of the ej. 

This situation prevails when ^ ej > 0. When ^ ej = 0, the driving parameter parameter is 

Normal forms can be used to compute the dimension of the preimage of a singular point. Fixing 
the values of the conserved quantities amount to fix the values of the quadratic terms. For an elliptic 
term we have 

p) + q) = hj 
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Figure 1: The solutions of eq, (|32l) when ej = — 1. 
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and the preimage is a circle. If hj = this circle degenerates to a point. In the case of a focus-focus 
singularity we have 

PjQj +Pj+iqj+i = hj 
I'j'U ■ ■ ~ Pj- ■<!;, = hj+i 

Solving in pi for instance we find 

_ hjqj + h j+ iq j+ i _ hjq j+1 - hj+iqj 

Pi ~ 2 , 2 ' Vj+\ ~ 2 , 2 
9,- + Qj + 9J+i 

so we have a two dimensional preimage. If hj = hj+i = 0, we have two planes qj = qj+i = or 
Pj = Pj+x = which intersect in one point. The preimage is a pinched torus. 



5 Monodromy. 

Let us go back to eq. ([3T1) . We have 

2 / \ \ , v- «( A ) 



Of course 



\Aj( x - e i) j a'(nj)(X - Hj) 

Uj(^-Pj) 



a(X) = 2 



n,(A 



We want to study the roots of the polynomial Q2n+2(A) when we perform a small circle around the 
singularity in the space Hi i.e. we want to examine the motion of the branch points of the spectral 
curve (to be defined in the next section, see eq. (|34p ) when we turn around a singularity. To zeroth 



order in the small deviations b,b,s- away from a critical point, the polynomials B{\) and C(A) 



vanish, whereas A(X) can be replaced by the fixed function a(A). Therefore: 

£?2n+2(A) 



n,(A-e 



j i 



\ - 



a\X) 



and Q2n+2(A) has double zeroes at the roots \Xj of the classical Bethe equation a(fij) = 0. Let us 
now pick a phase-space point close to the critical point. The corresponding invariant polynomial 
£?2n+2(A) is expressed in terms of the normal coordinates B(hj), C(fij) through eq. (f33|) Consider 
a root \ij . The branch point is moved at fij + Sfij . Inserting into eq. (j33[) we get the equation 
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so the leading variation of the root fj,j is due to the normal mode B(fj,j),C(fj,j) only. The other 
modes contribute in subdominant terms. This is Krichever's result |15j . 

If the root fij is real, we see that that it splits into a pair of complex conjugated roots, and 
the splitting is in the imaginary direction. To first order the deformation space is one dimensional, 
because, as we have seen, B(/j,j) = C(fij). 

If the root fj,j is complex then (5fij) 2 is complex. 

{5 H f = - (^) 2 B(H)C(H) = ^ {P fo^-focus +ip fo C us-focus ) 

The deformation space to leading order is now two dimensional. When we perform a small loop 
around the singularity in //-space, the branch points perform half a turn in opposite directions and 
at the end of the process they are exchanged. This is the basis of the interpretation by Michele 
Audin P2] of the monodromy [12\ [T7\ [T5| |2"U] using Picard-Lefschetz theory. 



6 Riemann surfaces and integrability. 

In this section we recall, in the example of the Jaynes-Cummings-Gaudin model, the algebro- 
geometric solution of classical integrable models [21, 22, 23, 24J. We first introduce the spectral 
curve. At each point of the spectral curve we can associate an eigenvector of the Lax matrix. When 
properly normalised the components of this eigenvector are meromorphic functions on the spectral 
curve. The poles of these meromorphic functions are coordinates on phase space. We compute their 
Poisson brackets. The image of the divisor of the poles by the Abel map is a point on the Jacobian. 
The motion of this point under the Hamiltonians of the system is linear. We deduce from this that 
when the spectral curve is non degenerate, the image of the moment map has maximal rank. 

The spectral curve is a curve in C 2 defined as det(L(A) — fj) = or: 

r : -A 2 {\) -B(X)C(X) = (34) 

Its importance arises from the fact that it is invariant under the time evolution, so it encodes 
knowledge of all the commuting integrals of motion. Equivalently, a spectral curve is attached to 
each fiber of the moment map. In the Jaynes-Cummings-Gaudin model, its equation reads: 

V 2 = 4A 2 + W n+1 + 2 £ A + £ 

We will consider the system reduced by the symmetry generated by H n+ \. It acts on the Lax 
matrix by conjugation by a diagonal matrix and obviously leaves the spectral curve invariant. 
Hence the dynamical Hamiltonians are Hj, j = l,---,n, and H n+ \ should be considered as a 
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constant parameter. Note that the dynamical Hamiltonians Hj appear linearly in the equation of 
the spectral curve: 

9 

r : R(\,n) = Ro(\,ii) + J2R j {X)H j = (35) 

3=1 

where: 



RjW = 1 . ^o(A,m) = -fJ 2 + 4A 2 + 4:H n+ i + — 

A - €j . (A - €j) 

As we already noticed 



A 2 (X) + B(X)C(X)- Q ^+^ 



where Q2n+2(X) is a polynomial of degree 2n + 2. Defining y = nYlj(X — tj), the equation of the 
spectral curve becomes 

V 2 = Q2n+2{X) 

which is an hyperelliptic curve of genus g = n. The dimension of the phase space of the model is 
2(n + 1). However, we reduced it by the action of the group of conjugation by diagonal matrices 
which is of dimension 1. Hence we confirm that 

g = ^dim A4 rcduccd = n 



At each point of the spectral curve, we can solve the equation 

(L(X) - a*)* = o 

Normalizing the second component (instead of the first one, for later convenience) of * to be 1, we 
find: 



1 J ' ri C(X) 

showing that tpi is a meromorphic functions on the spectral curve. The poles of VP at finite distance 
are located above the zeroes of C(A). Note that if C(Afc) = 0, then the points on T above have 
coordinates fi^ = ±A(Xk)- The pole of VP is at the point fi^ = A(Xk), since at the other point, 
= — A(Xk), the numerator of V'l has a zero. 



Recalling that 



™ = *+EiV*fe&^ <36) 

we see that indeed the eigenvector has n = g dynamical poles. 
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At infinity, we have two points: 

Q± : 11 = ±2\{1 + 0(\- 2 )) 

Remembering that 

A(\) = 2\ + 0(\- 1 ), C(\) = 2b + 0(\- 1 ) 
we find the following behavior of the function ip\ at the two points Q± 

Q+ :Vi = \ X + OiX- 1 ), Q_ : Vi = 0(A _1 ) 

showing that the eigenvector has a pole at Q + and a zero at Q_ in agreement with the general 
theory. By the Riemann-Roch theorem, the meromorphic function ipi exists and is unique. 

We can reconstruct the phase-space point of the reduced model (where fln+i is a fixed param- 
eter) from the knowledge of the 2n variables (Xk,fJ>k)- From eq. ([36]) . we get: 

S t = 2b Iff=ifo- A *) (37) 

\\i.-j u .i f ' ! 

The components can be obtained from the n equations: 



— 2Xk = 



after inversion of a Cauchy matrix. This analysis shows that the Lax matrix can be reconstructed 
once we know the coordinates (Afc,//^) of the poles of the eigenvectors. Hence (Xk,^k) can be 
considered as coordinates on the (reduced) phase space. We now compute the symplectic form in 
these coordinates. 

From the constraint (s|) 2 + s j^ s J = s 2 we can eliminate the variables sj . Remembering the 
Poisson bracket {sj, s^} = isj, we can write the symplectic form as 



n = -i5b A 5b + i — ^- A 6s j 



From eq. (|37j) . we get: 



therefore: 



6s t _ Sb sr 6Xk 



at b ^ X k 



= A 5b + i j- A 2^ <5s* + z 2^ 2^ 



Afe 
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But . . 

SA(X k ) = 25X k + V - Sj 5\ k 

j *k~ tj (Afc - ej) z 

and 



n = i^- a 

b 



(5(66) + Ssj] +iJ2 6X * A 6A ( x k) 



Finally: 

Q = iS log 6 A 5H n+ \ + i (JA^ A 5fi k 

k 

This shows that the variables (X k ,fJ, k ) are canonically conjugate. The above calculation is valid 
before the symplectic reduction by H n+ \. This expression for the symplectic form confirms that 
the separated variables are invariant under the diagonal group action 

{H n+l , X k } = 0, {H n+1 , iM k } = (38) 

so that they can be used as coordinates on the reduced phase space. 

We have found that T is of genus g = n and there are exactly g commuting Hamiltonians Hj. 
Moreover ^ has exactly g dynamical poles. The curve is completely determined by requiring that 
it passes through the g points (Aj,//j), i = 1, ■ ■ ■ ,g. Indeed, the Hamiltonians Hj are determined 
by solving the linear system 

E A^b"^ = \» l " 2X1 " 2Hn+1 (X k = e,) 2 (39) 

whose solution is 

H = B~ l V (40) 

Here the matrix B k j is the Cauchy matrix 

Bki = (41) 

*k — £j 

and V k is the right hand side of eq. ([39|) . 

We now compute the equations of motion of the A^'s. One has 

d ti X k = {H h X k } = -Y J {B u 1 V h X k } = -Y,B u 1 {V u Xk} = -B7 k 1 {V k ,X k } 
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where in the last equality, we used the separated structure of the matrix B and the vector V to 
suppress the sum over I. Explicitely 

dti^k = ~i{B~ 1 )ikfJ'k ( no summation over k) (42) 



In order to write the equations of motion eq. (|42p we need to invert the Cauchy matrix. We find 

( )jp "n^(^-^)n^(A P -Ao (43j 

Hence 



~ x , VQ2n+2(Afc) H#fc( e i ~ Aj 

n^ fc (Afc-A/) n#ite-ej) 

where the — i comes from the symplectic form. 
Note that we can also write equivalently 

^-Bjfc — S^Afc = -iSij 
k ^ k 

where 

/\ \ B kj __ n^-(Afc-Q) 

^'(AfcJ - - — ====■ 

^fc V < 32n+2(Afc) 

but aj(\)dA are precisely a basis of holomorphic differentials. Hence we have 

^2 dt ^ k a i( X k) = ( 45 ) 



Define the angles as the images of the divisor (A/j,/^) by the Abel map: 

Qj = }2 <Tj(\,(Jt)d\ 
k •* 

where <7j(A, fi)d\ is any basis of holomorphic differentials. This maps the dynamical divisor 
(Afc,/Ufc), k = 1" -n to a point on the Jacobian of V. We can now prove the fundamental theo- 
rem of classical integrable sytems 

Theorem 1 Under the above map, the flows generated by the Hamiltonians Hi are linear on the 
Jacobian. 
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Proof. We want to show that the velocities d^Oj are constants, or 



k 

but this is eq. ([4"5j) . 



ste 



When the divisor (X k ,fM k ) is in general position, the Jacobi inversion theorem [25J implies that 
the matrix crj(Afe) is invertible and therefore so is the matrix c^Afc. An important consequence of 
this fact is that the projections on the reduced system of the flows t{ are all independent as long as 
the spectral curve is non degenerate. In that case, the moment map is of rank (n + 1) if the orbit 
of H n+ \ is one dimensional, or of rank n if the orbit of H n+ \ is of dimension zero. But the flow 
generated by H n+ \ is just a phase 

b(t n+1 ) = e^b(Q), sf(t n+1 ) = e ±u ^sf(0) 

For the orbit of H n+ i to be of dimension zero we must have 

6 = 6 = 0, = 0, s j = se ji e j = ±1 

and these are precisely the critical points where the rank of the moment map is zero. Outside these 
points the rank of the moment map of the full system is equal to the rank of the moment map of 
the reduced system plus one. 

Let us assume that Q 2 n+2(A) has a double real root at A = E. 

Q 2 n+ 2 (A) = (A - E) 2 Q 2n (X) 

This means that ^4 2 (A) + I?(A)C(A) has a double root. But for real A one has C(A) = B*(X) so 
that A(X), B(X), C(A) must all vanish at A = E. 

Recalling eq. (|36j) . this meang^l that one of the separated variables, say X\(t), is frozen to E. 
This is compatible with eq. (!44p which becomes 



8 u X k = -i v TT-frT \~( e i ~ E )i k = 2---n (46) 

These flows are not independent however. Because of the identities 

n p 

^Zrr-r 1 v = ' o< P <n-2 



2 If b(t) = b(t) — 0, then B(\) or C(A) are identicaly zero. The equation of motion dtjb = igsj implies s^(t) = 
and therefore Sj(t) — sej. Hence this corresponds to the singular points. 



17 



we have the relation 



1 



d ti X k = 0, k = 2---n 



-E 



On this submanifold the rank of the moment map of the reduced system is therefore n — 1, and 
the moment map of the full system has rank n. By requiring more and more double zeroes and 
freezing more and more variables we construct the different strata of the moment map. This 
connection between degeneracies of the spectral curve and singularities of the momentum map has 
been exploited before for several integrable systems such as spinning tops [8]. 

In the case of a double complex root E, we do have solutions where one A& is frozen to E, but 
these do not exhaust all the corresponding fiber of the moment map. A detailed example of this is 
given in section [7] below, in the case of a system with one spin (n = 1). 

We arrive at the interesting conclusion that the degeneracies of the moment map and the 
degeneracies of the spectral curve are intimately related. Recall that the spectral curve reads 



The allowed real polynomials Q2n+2(A) of degree 2n+2 are characterized by the following conditions 
in the above expansion: 



These are n + 2 conditions on the 2n + 3 real coefficients of Q2n+2(A). The remaining n + 1 
coefficients are precisely the n + 1 Hamiltonians Hi . The degeneracies we look at are of the form 



where the coefficients aj , bj are real. The integer < r < n + 1 will be the rank of the moment 
map. To make the decomposition unique we can always impose a ra +i_ r = 1 so that we have 
(n + 1 — r) + 2r + l = n + r + 2 coefficients on which we impose n + 2 constraints. Hence the leaf 
of rank r is of dimension r. We remark that the conditions we have to impose appear as linear 
equations on the bj so that we always start by solving them. Another consequence of this remark 
is that if 2r > n + 1 it remains 2r — n — 1 free coefficients bj which enter the problem linearly. 




(47) 



The coefficient of A 2 is equal to four, 
The coefficient of A is equal to zero, 
The coefficients of the double poles are s 2 
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Hence the strata of rank r contain linear varieties of dimension 2r — n — 1. In the cases of rank 
and rank 1, the genus of the spectral curve is zero. 

Strictly speaking, the rank can vary along a given fiber of the moment map. For example, let 
us consider a focus-focus critical point in an integrable system with two degrees of freedom. As 
we have already mentioned in section [4] the preimage of the critical value of the moment map is a 
torus which is pinched at the critical point. Therefore, the rank of the moment map for a generic 
point on this pinched torus is two, but it falls to zero at the critical point. The above prescription 
of freezing n + 1 — r separated variables A& on double zeroes of Q2n+2 picks configurations which 
minimize the rank on the fiber of the moment map defined by the spectral curve. 



7 The one-spin model. 

We now study the example of the one-spin system from the point of view of the degeneracies of 
the spectral curve. This model is very well known in the physical literature [3] but also appeared 
recently in the mathematical literature [13] . In the one-spin case, the Hamiltonians read 



Hi = 2eisf + bsf + bs± 
H 2 = bb + sf 



Recall again the spectral curve eq. H19H which reads in this case : 

^% = 4A 2 + 4F 2 + 2-*- + (48) 
{X-eiY A-ei (A-ei) 2 

We want to see when the spectral curve degenerates. 
7.1 Rank 0: 

As we have seen the singular points are given by b = b = 0, sf = 0. Hence we have two points 

s\ = es, e = ±1 
The corresponding values P = (Hi, H 2 ) are 

Pi(t) = [2e lS ,s] 
P 2 (i) = [-2e lS ,-s] 

We can recover this result by analysing the spectral curve. The most degenerate case is when 
Q4(A) is a perfect square. Assuming 



Qi(\) = (a 2 A 2 + aiA + a ) 2 , a 2 / 
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We have three coefficients Oj on which we impose three conditions hence they are completely 
determined. We find 

a 2 A 2 + aiA + a = 2 ^A 2 - eiA + ^ejj 

Comparing the partial fraction decomposition of both sides of eq. (|48p . we find the corresponding 
values of (-f/i, P^)- They are precisely Pi and P%. Hence the points of rank zero are the only points 
where the spectral curve is totally degenerate. 

To determine the type of the singularities, we look at the classical Bethe equations which read 

se o 

2(i H = -O- 2/T - 2ei(i + se = 

(i - ei 

The discriminant of this equation is e 2 — 2se. So, when the spin is down (point P 2 ), the discriminant 
is positive, the two classical Bethe roots are real and this is an elliptic singularity, in agreement 
with the general analysis of section HJ When the spin is up (point Pi) we have real roots when 
e\ > 2s (i.e. the singularity is elliptic in this case), and a pair of complex conjugate roots E,E 
when e 2 < 2s (i.e. the singularity is focus-focus in that case). 



-1.4 1.4 €\ 
■ ■ ■ > 

Stable Unstable Stable 



Figure 3: The domains of the coupling constant ei, (s = 1). 

Let us now discuss the fibers of the moment map over the critical values Pi and P 2 . In the stable 
case, such fiber is reduced to the critical point. But in the unstable case (focus-focus singularity), 
this fiber is a two-dimensional torus pinched at the critical point. After the symplectic reduction 
associated to P 2 , the pinched torus is conveniently described as a finite arc in the complex plane 
of the separated variable Ai. The spectral polynomial can be written as: 

Q 4 (A) = 4(A-P) 2 (A-P) 2 (49) 

where: 

ei = E + E, EE= S - 
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The separated variable Ai is defined by: 



C(A) = 26- = 26 + 



A — ei A — ei 

so sf is expressed as: 

s + = 26(ei - Ai) (50) 

The conjugated variable n± is equal to 

/ii = A(A X ) = 2Ai + 



A-ei 
so that: 

sf = (ei-Ai)(2Ai- M i) (51) 
Since (Ai,/ii) belongs to the spectral curve, we have: 

,(A 1 -S)(A 1 - J E) 



Mi = ±2- 



(Ai - ei) 



Choosing the plus sign leads to sf = s which is the unstable point. So we now choose the minus 
sign, which gives: 

s\ = -2 (Ai(Ai -E-E) + (\ 1 - E){Xi - E)) (52) 
The equation of motion for the flow generated by H\ is: 

d tl X 1 = 2i{\ l -E)(X 1 -E) (53) 

whose solution is: 

E-EX X = A 2i(E-E) tl 

1-X 

From the expression of sf eq. (j50|) and the equation of motion 

d tl b = -is+ (54) 

we deduce: 

' ' ~ 1 - X ' 1 ~ (1 - X) 2 

From eq. ([52|) we get: 

sf = 2EE-A(E-E) 2 * 

Now sf should be real, which is equivalent to (X - X)(XX - 1) = 0. With X = A x e 2i ^ E ~ E)tl it is 
impossible to have XX = 1 for all t\. Then X = X , which imposes A\ = A\ so that A\ is real. Its 
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absolute value can be absorbed in the origin of time t±. Only its sign matters. To determine this 
sign, let us impose the constraint that the length of the spin is constant. This gives the condition: 

b b = 4A 1 (E-E) 2 

Then: 

( S *)2 + s + s - = AE *e 2 = s 2 , 

We see that A\ should be negative. As a consequence, Ai runs along the line interval joining E and 
E. Here, we see explicitely that freezing Ai on E (or on E) gives only the critical point (t) which 
is only a tiny part of the fiber above the critical value Pi in the focus-focus case. Note that on this 
fiber the rank of the moment map is two, excepted fot the critical point (f) where it vanishes. 



7.2 Rank 1 

We assume next that 

Q 4 (A) = (A + a ) 2 (M 2 + &iA + & ), b ? + ( 55 ) 
where ao and bi are real. We denote: 

ao = -x — ei 

Imposing the three conditions on Q4(A), we can determine bo,bi,b 2 in term of x by solving linear 
equations. We find: 

b 2 = 4, bi = -4x, b = 4 % 



H 2 = - 



x 4 - 2eix 3 - 4s 2 
2x 

3x 4 -8eix 3 +4e 2 x 2 -4s 2 



4x 2 

and we see that: 

OHi = x dlh (56) 
dx dx 

Since we are on a rank one line, the derivatives of the functions H\ , H 2 with respect to any coor- 
dinates X and Y on phase space, evaluated on the line are proportional: 

d ( HA d ( HA dxHx dyH x 

oc — — TT ==>- — — = — — — = X 



dX\H 2 ) dY\H 2 J d x H 2 d Y H 2 

In particular we have: 

9Hi 8H 2 + 
-db =X ^b^ S ^ =xb 
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Inserting this relation into the definitions of H\ and H2 gives: H2 = bb + sf and H\ = 2e\s\ + 
2xbb. Therefore: 

sf = ~-x(x - 2ei) 

l h _ ( 2s -x 2 + 2e 1 x){2s + x 2 - 2e±x) 
Ax 2 

Notice that when 

x 2 - 2e x x + 2es = 0, e = ±1 (57) 

we have sf = es and bb = and this corresponds to the points Pi and P%. In fact there is a simple 
relation between rank and rank 1. Rank 1 spectral curves degenerate when the polynomial 
62 A 2 + b\ A + 60 has a double real root. Its discriminant is: 

,2 „, (x 2 - 2eiX - 2s)(x 2 - 2eiX + 2s) 

b 2 - 4b b 2 = 16 ^ ^ - = -64 bb < 

It vanishes precisely when eqs. ([5T|) are satisfied. These equations are nothing but the classical 
Bethe equations eq. ([26|) as can be seen by setting 

x = 2fi 

Remark that the boundary of the image of the moment map is obtained for x real, because Hi 
and H2 are both real. Hence the elliptic points are on the boundary, but the focus-focus points 
correspond to complex x and therefore are not on this boundary (see fig. [7] below). To end the 
characterization of the boundary of the moment map we have to determine the range of x. The 
physical constraints are — s < sf < s and bb > 0. These two conditions are both equivalent to: 

-2s + e\ < (x-ei) 2 < 2s + e\ (58) 

In this inequality we recognize the discriminants e\ ± 2s of the classical Bethe equations, so that 
we have to distinguish between the stable and unstable case. 



7.2.1 Stable case. 

When e 2 > 2s the two sides of the inequality are positive and put bounds on x. The two equations 
determining the boundary values of x are precisely eqs. (|57p . and we are in the case where all their 
roots are real. The allowed range of x is as in fig.Q. From this we can construct the image of the 
moment map. It is shown in fig.© with the edges labelled according to the ranges of x. The point 
P2 (spin down) is the green point and the point Pi (spin up) is the cyan point. They are both on 
the boundary of the image. 

The swallow tail is not part of the image of the moment map. It is composed of the points 
{H\{x) 1 H2{x)) when x runs in the interval [ei — \J — 2s + e'f, e\ + \J — 2s + e'f\. The two cusps 
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swallow tail 


— • • — 

2 3 


> 



Figure 4: The range of x. The allowed segments are labelled 1, 2, 3. The points of the same color 
are mapped to the same singular point in the bifurcation diagram. Notice that x = is mapped to 
infinity, (ei = —2). 
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/ 3 




1 










6! = -2.0 







Figure 5: Image of the moment map in the stable case. The swallow tail is not part of the image 
of the moment map. (ei = —2). 



correspond to values x c of x such that dH\jdx and dH2/dx both vanish. Let us write x = x c + u. 
Then, taking into account eq. (|56|) . we have the following Taylor expansions for small u: 

5H 2 = ^u 2 + ^u 3 + 0(u 4 ) (59) 
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5Hi = x c 5H 2 + -u 3 + 0(u 4 ) (60) 
where 5 Hi stands for the the small variation Hi(x c + u) — Hi(x c ). From eq. (|59|) . we have: 

/ 2 \ 1/2 
u = ±(-\8H 2 \\ +0(5H 2 ) 



Therefore: 



/ o \ 3/2 



8H X = x c 6H 2 ± - l-\SH 2 \ ) + 0{{8H 2 f) (61) 



which gives the leading shape of the cusp near x c . We see that the two branches are found on 
opposite sides of their common tangent, as shown on fig. [5j 



7.2.2 Unstable case. 

When e\ < 2s the left side of the inequality ([58]) is always satisfied. Among the equations eqs. ([57]) 
determining the bounds of x, one has two real roots and the other one has two complex conjugate 
roots. The range of x is as in fig.©. The image of the moment map is shown in fig. ([7]). The point 
P 2 (spin down) is the green point and is a vertex on the boundary. The point Pi (spin up) is the 
red point. It is in the interior of the image of the moment map. 








X 


1 


2 


> 



Figure 6: The range of x. The green points correspond to the stable point P 2 and are solutions of 
the equation x 2 — 2t\x — 2s = 0. Since points of the same color are mapped to the same point by 
the moment map, and zero is mapped to oo, the boundaries of the image of this map are easy to 
reconstruct. {e\ = —0.707). 
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Figure 7: The image of the moment map in the case of one spin, with one unstable critical point. 
The green point is the stable (elliptic) point P 2 . The red point is the unstable (focus-focus) point 
Pi. It is in the interior of the image of the moment map. (ei = —0.707). 



8 The two-spins model. 



For the two spins model, let us write explicitly the Hamiltonians: 

Hi = 2eisf + bsf + bs^ + Sl ' S2 

ei - £2 
si ■ s 2 



where 



H 2 = 2e 2 s z 2 + bsj + bs 2 

ei — £2 

H 3 = bb + sf + s^ 



si ■ s 2 = s\ s\ + -(s 1 4 + sf s 2 ) 



The spectral curve reads in this case: 



<2e(A) 2H i 2H 2 s 2 s 2 

~ 4A + 4 #3 + - + i — + - — + - — (62) 



(A-ei) 2 (A-e 2 ) 2 A - ei X - e 2 (A-ei) 2 (A - e 2 ) 

Let us now use the degeneracies of zeroes of Qq(\) to study the rank of the moment map. 
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8.1 Rank 

Again, the singular points are given by b = b = 0, sf = sf = so that we have four critical points: 

sf = ±s, si = ±s 
The corresponding values P = (Hi, Hi, H-$) are: 



A(t,t) 

P2(ti) 

P 4 (U) 



S 2 S 2 

h2eis, 2e 2 s . 

ei - e 2 ei - e 2 

„2 „2 



2s 



2eis 



ei - e 2 ei - e 2 



-2ei s 



ei - e 2 



ei - e 2 ei - e 2 



2e 2 s, 



+ 2e 2 s, 



2eis, -2e 2 s 



-2s 



ei - e 2 



(63) 
(64) 
(65) 
(66) 



We can recover this result by analysing the degeneracies of the spectral curve. The most 
degenerate case corresponds to a polynomial Qe(X) of the form 

Qe(A) = (a 3 A 3 + a 2 \ 2 + oiA + a ) 2 

The four conditions we have to impose on Qq(X) determine completely the four coefficients ctj. We 
find: 



a 3 A 3 + a 2 A 2 + oiA + a = 2 (A 3 - (e x + e 2 )A 2 + (eie 2 + -(e a + e 2 )l A - -(e a e 2 + e 2 ei 



(67) 



where ej = ±1. The values of the energies are exactly eqs. (j63H66j) . 

In order to determine the type of the singularities we write the classical Bethe equations: 



sei se 2 

2/i H 1 

M - ei /x - e 2 







These are polynomial equations of degree 3. The number of real roots is determined by the sign of 
the discriminant: 

Disc(ei, e 2 ) = -4e 2 e^ (e 2 - eif + 4 (e 2 - ei) (2 + e 2 ei - 2eieie^ + 2e 2 e 2 e 2 - eie 2 e 2 - 2e 2 ef) s 
- (e| + 20eie 2 e| - 8e| + 8eie 2 - 38eie 2 eie 2 + 8e x e 2 - 8e 2 + 20eie 2 e 2 + e 2 ) s 2 
+2(e 2 + ei) 3 s 3 

If Disc(ei, e 2 ) < we have three real roots (elliptic point) and if Disc > we have one real root 
(focus- focus singularity). We restrict ourselves to the case where ei and e 2 are both negative. We 
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Figure 8: The domains indicating the signs of the various discriminants, for three critical points 
of the system with two spins. Curves correspond to the vanishing of one discriminant. The color 
of these curves and of signs inside the domains label the critical points according to the following 
code: red stands for up-up, green for up-down, and blue for down-up. (s = 1). 



easily see that Disc(ei = — l,e2 = —1) is always negative, hence the configuration with two spins 
down is always stable. The signs of the other three discriminants is shown in fig. dHJ) - 
In region V, for instance, we have two unstable points and two stable points. 

(t, t) : one rea l ro °t (unstable) 

(t, 4) : three real roots (stable) 

(4, t) : one real root (unstable) 

(4, J.) : three real roots (always stable) 
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8.2 Rank 1 

We now set 



Qe(A) = (A 2 + aiA + a ) 2 (b 2 \ 2 + h\ + b ) (69) 



We have five coefficients and four conditions on Qq. Hence we have a dimension one manifold 
of solutions. The coefficients bj are completely determined and there is one constraint between 
(ao,a\). Note that the case of rank is obtained as a special case of rank 1, when the polynomial 
6 2 A 2 + bi\ + 60 has a doubly degenerate root. Let us parametrize (ao,ai) in terms of (x,y) as 
follows: 

ei V + £2 x - 2ei e 2 y + x-2e 2 -2ei 
ao = g , ax- - 

Imposing the coefficient of A 2 we find: 

62 = 4 

Imposing the vanishing of the coefficient of A, we find: 

61 = -4 (y + x) 

Imposing that the coefficient of the double pole at A = e 2 is s 2 we find: 

4 (e 2 y 3 + e 2 xy 2 -e 2 2 y 2 + s 2 ) 
bo — 5 

yZ 

Imposing next that the coefficient of the double pole at A = e\ is s 2 we find the constraint: 

Si: s 2 (y 2 -x 2 ) + (e 1 -e 2 )(x + y-e 1 -e 2 )x 2 y 2 = (70) 
Then H U H 2 , H 3 are given by eqs. (|71|72)73p . 

Hi = — l -^{{e 2 -e l ){-2xy 4 -2,x 2 y 3 -2 (-3e 2 + ei) xy 3 - x 3 y 2 ) 

2 (e 2 - ei) y l 

-(e 2 -e 1 ) (-4e 2 + 2e 1 ) x 2 y 2 + 4 (e 2 - e x ) 2 (-e 2 - d) xy 2 -2s 2 x?/ 

+4 (e 2 -ei) s 2 x) (71) 



„ -(e 2 -ei) y 4 -(e 2 - ei ) xy 3 -(e 2 -ei) (-2e 2 ) y 3 + 2s 2 x + 4e 2 s 2 -4eis 2 

2 (e 2 - ei) y 



-3 j/ 4 - 6 x y 3 - 4 (-2 eg) - 3 x 2 y 2 - 4 (-e 2 - eQ x t/ 2 - (-2 e 2 ) 2 7/ 2 + 4 s 2 

^ 3 = ly~2 ^ 
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These are the parametric equations of the lines of rank 1. The parameters x and y are tied together 
by relation eq. ([70|) . 

Let us now define the total derivative with respect to x by: 



Then we can compute: 



d d dy d dy d x S\ 

dx dx dx dy' dx d y Si 



^h_ x <^h _ V + 2(ei - e 2 ) dy 

dx dx (e\ — e 2 )xy 3 dx 

dH 2 dH 3 2 dy s 

dx dx (e% — e 2 )xy 2 dx 



hence, on Si, we have: 



dHi dH 3 dH 2 dH : 



x -J - ) ~j = V 



3 



dx dx dx dx 

Because we are on a rank one line this relation is true for any derivative on the line. Considering 
the derivative with respect to b we find 

Si = xb, S2 = yb 
and considering the derivative with respect to sf we get: 

1 1 

2sf H {-xs z 2 + ysf) = -x 2 + 2eix, (~xs 2 + ys\) = xy 

e\ — e 2 ei — e 2 



from which we deduce: 



~(x + y-2e 1 ), s z 2 = - y -{x + y-2e 2 ) 



and: 

- = {y 2 + xy- 2e 2 y - 2s)(y 2 + xy - 2e 2 y + 2s) 

Note that the discriminants of the two second degree polynomials which appear in the factor- 
ization (|69p of Qq(X) are: 



,2 Ah , 16 (y (x + y)- 2e 2 y - 2s)(y(x + y - 2e 2 y + 2s) 

\ - 4b b 2 = = -t 

aj-Aao = \((x + y) 2 -4(e 1 -e 2 )(x-y)+4(e 1 -e 2 ) 2 ) =]a 
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The first discriminant vanishes when Qq(X) has three double roots, that is when the rank drops to 
zero. In this case, inserting the parametrization 



x = — s 



y = -s 



V - e 2 



into eq. (|70p . we see that it turns into the classical Bethe equation for /i. The variables x, y being 
real by definition, the real solutions of the classical Bethe equations correspond to points on the 
curve Si. 

When A > 0, Qe(X) has two real double roots, which become a pair of complex conjugated 
double roots when A < 0. Most likely, the sign of A governs the nature of the preimage of points 
lying on such a line of rank 1. To check this, one would need to compute the corresponding normal 
forms, which we have not done yet, the normal forms discussed previously being specialized to 
the special case of rank zero singular points. We leave this generalization as a subject for future 
work. We expect that points where A = correspond to qualitative changes in the topology of the 
pre- image. They are determined by the following factorization for Qe(X): 



Q 6 (X) = A(X-a ) 4 (X 2 + b 1 X + b ) 



As before, we express the known constraints on this polynomial, which gives: 



4 (-a + ei) 4 (£> + h ei + ef) - (ei 

4 (-a + e 2 ) 4 (6o + h e 2 + e|) - (ei 
-16 ao + 4&i + 8ei +8e 2 







6 2 ) 2 * 2 = 



6 2 ) 2 * 2 = 



Solving for 6q and b\ in the first two equations gives: 



fro = ei £2 



(-ef + eie 2 ) s 2 (-ei e 2 + e§) s 2 
4 (o - e 2 ) 4 4 (o Q - ei) 4 



and 




(-ei + e 2 ) s 2 (-ei + e 2 ) s 2 
4 (a - ei) 4 4 (a - e 2 ) 4 



Plugging these values in the third equation determines qq through: 




= 
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8.3 Rank 2 

We now set: 



Q 6 (A) = (A + a ) 2 (6 4 A 4 + 6 3 A 3 + 6 2 A 2 + b ± X + 6 ) (74) 



The four constraints are linear equations determining 64,63,62,^1- It remains two parameters 
ao,6o- The rank 2 manifolds are two-dimensional surfaces. Remark that 60 appears linearly. The 
four conditions on Qq read: 

64 = 4 

(ei + a ) 2 (64 e\ + 63 e\ + b 2 e\ + 61 e x + b ) = (ei - e 2 ) 2 s 2 
(e 2 + a ) 2 (64 e 2 + 63 e 2 + 6 2 e 2 + 61 e 2 + 6 ) = (ei - e 2 ) 2 s 2 
2 64 e 2 + 2 6 4 ei + 2 a 64 + 63 = 

We can solve them in terms of ao and 60: 

6 4 = 4 

63 = -8 (e 2 + ei + oq) 

, (£2-ei)s 2 (e2-ei)s 2 4 e x e% + a (8 e : + 8 e 2 e 2 ) + 12 ef e 2 , + 4 ef e 2 + b 
"2 — : -o ; 79 + 
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e 2 (e 2 + a ) ei (ei + a ) e i e 2 

(e| ~ d e 2 ) s 2 _ (ei e 2 - e 2 ) s 2 _ 4e 2 e| + 4ef e| + 8a e| e| + 6 e 2 + 6 ei 
ei (ei + a ) 2 e 2 (e 2 + a ) 2 £ i e 2 



Once these constraints are implemented, we can write the Hamiltonians on the faces of the image 
of the moment map, which are therefore parametrized by ao and 6 : 



Hi = 



(e 2 2 + ej) s 2 + a (8 ef e\ + 8 e\ e 2 - 2 6 e x ) + a 2 , (4 e 2 e 2 + 16 ef e 2 - 6 ) + 4 e 4 e 2 , + 8 a 3 e\ e 2 - 6 e 2 



2ei e 2 — 2e 2 e 2 



( 1 s 2 (' 1 '2 - f ■') s 2 -s 2 



+ - ; ^ + 



Ho 



e 2 (e 2 + a ) 2 e 2 (e 2 + a ) 2 e i + ao 

( e 2 + e i) g 2 + Qo (8eie| + 8e 2 e^ - 26 g 2 ) +4e 2 e 2 t + a 2 > (l6ei e\ + 4ef e 2 . - 6 ) + 8ag ei e| - bp e| 

2 ei e| - 2 e? e 2 

s 2 (e 2 . - ei e 2 ) s 2 e 2 s 2 



e 2 + a 2ei(ei + a ) 2 ei (ei + a ) 

(e 2 -ei)s 2 (e 2 -ei)s 2 a (8ei e 2 , + 8ef e 2 ) + 4ef e\ + 12a 2 , e x e 2 - 6 



4e 2 (e 2 + a ) 4ei(ei + a ) 4eie 2 
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Remark that bo enters the formulae linearly so that the rank two faces are ruled surfaces. To find the 
intersection between two faces, we go back to eq. (|69p . As we have seen, the discriminant b\ — 46o&2 = 
—64 bb is zero at a critical point and negative as soon as we leave this critical point. So the double 
real root at a critical point splits into a pair of complex conjugate roots as soon as we move along 
the rank 1 line. If the other discriminant a\ — 4ao = (1/4)A is positive, the polynomial A 2 + a±X + ao 
has real roots ot\ and «2 5 and we can recast Qe as Qg(X) = (A — ai) 2 ((A — a2) 2 (^2A 2 + b\\ + bo)). 
Expanding the second factor (of degree four), we obtain an expression of the form (|74|) with real 
coefficients, which satisfies all the constraints. This establishes that the line of rank one, defined 
by eq. (|69p is included in a face of rank two, as long as the roots a± and «2 are real, that is when 
A > 0. When A crosses zero to become negative, the line of rank one starts to leave the face and 
goes inside the image of the moment map. The corresponding roots form a complex conjugated 
pair of double roots. 

We have also established that the points of rank zero are on the faces as well as the edges of 
rank one, provided A > 0. This last condition excludes the lines of rank 1 which go in the interior 
of the image of the moment map when A < 0. 
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Figure 9. The curve SI for a point in region /. The red and blue points correspond to the unstable points Pi 
and P3 respectively. The green and cyan points correspond to the stable points P2 and P4 respectively. Since 
points of the same color are mapped to the same point by the moment map, the edges of the image of the 
moment map are easy to reconstruct. The black curve corresponds to A = 0. The origin has A > 0. We see 
that the edge joining the unstable points is not on the faces, (ei = —2, £2 = —0.2). 




Figure 10. The image of the moment map (IMM) and the bifurcation diagram for a point in region /. It has 
four faces. The stable points are vertices of the IMM. The unstable points are on and inside faces. A line of 
rank 1 joins the red point to the blue point passing inside the IMM. (ei = —2, £2 = —0.2). 
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Figure 11. The curve SI for a point in region II. The red, green, and blue points correspond to the unstable 
points Pi, P2 and P3 respectively. The cyan points correspond to the stable point P4. Since points of the same 
color are mapped to the same point by the moment map, the edges of the IMM are easy to reconstruct. The 
black curve corresponds to A = 0. The origin has A > 0. We see that the edge joining the unstable points is 
not on the faces. Note that the arc joining P2 (green point) to the origin (point at infinity in the IMM) crosses 
the A = curve, (ej = -1.5, e 2 = -0.3). 




Figure 12. The IMM and the bifurcation diagram for a point in region II. The stable point is the vertex of 
the IMM. The unstable points are on and inside faces. A line of rank 1 joins the red point to the blue point 
passing inside the IMM. Another line joins the unstable point P2 (green) to infinity. Starting from P2, it runs 
inside the IMM and later becomes an edge between two faces as the discriminant A becomes positive. This is 
a consequence of the crossing between the A = curve and line 2 shown on the previous figure, (ei = —1.5, £2 = 
-0.3). 
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Figure 13. The curve Si for a point in region V. The red and blue points correspond to the unstable points 
Pi and P3 respectively. The green and cyan points correspond to the stable points P 2 and P4 respectively. 
Since points of the same color are mapped to the same point by the moment map, the edges of the IMM are 
easy to reconstruct. The black curve corresponds to A — 0. The origin has A > 0. We see that the two 
edges starting from the unstable points are not on the faces, but they will join the faces when A becomes > 0. 
(ei = -1.2, e 2 = -1.735). 




Figure 14. The IMM and the bifurcation diagram for a point in region V. The stable points are vertices of 
the IMM. The unstable points are on and inside faces. A line of rank 1 joins the red point Pi to the green 
vertex P 2 . In the vicinity of Pi, it runs inside the IMM. But as shown on the previous figure, A changes sign 
along this line (where the color changes from green to red) before reaching P 2 . When A is positive, this line 
becomes an edge separating two faces. Another line of rank 1 joins the blue point P3 to infinity, passing inside 
the IMM in the vicinity of P3, and becoming later an edge between two faces, (ei = —1.2, e 2 — —1.735). 
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9 Conclusion 



In this article we have analyzed the bifurcation diagram of the Jaynes-Cummings model. The 
use of Lax pair techniques proved to be very useful. The classical analogue of algebraic Bethe 
Ansatz allows a very easy construction of the normal forms near critical points. We have shown 
in this way that this model possesses singularities of the elliptic and focus-focus type only. This 
is an approach alternative to the one due to Krichever [15] and based on the spectral curve. The 
spectral curve however was a very powerful tool to draw the full bifurcation diagram as advocated 
by Michele Audin [8J. In the one spin case we get results in agreement with general considerations 
([13, 18J), while in the two spins case it exhibits quite a rich structure which calls for more detailed 
investigations. 

Along the open questions is the determination of normal forms along the lines of rank one, and 
the explicit construction of real solutions of the equations of motion along these lines. Another 
fascinating subject is the emergence of this classical geometry from the quantum system and Bethe 
equations [26J. We hope to return to these questions in future publications. 
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